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ABSTRACT 

7 

Spectral peaks corresponding to tfie spheroidal free periods of oscillation of the earth exist 
in the spectra of eight seismograms written at stations in different parts of the world shortly 
after the great Chilean earthquake of 22 May I960. These data, have been combined with those 
previously reported by various authors to obtain a very precise phase velocity vs period curve 
for Rayleigh waves in the period range of 200 to 3200 seconds. The observed spectral amplitudes 
lend some support to the assumption of amoving source, but they also indicate t hat the present 
theory is not adequate. The vertical motion is found to be symmetric with respect to reflections 
through the pole. 


Introduction 


* 


Ai7fy-o< 


The seismograms written by vertical component seismographs at eight stations 
around the world shortly after the great Chilean earthquake of May 22, 1960, 
have been Fourier analyzed. The resultant spectra show peaks corresponding to the 
free oscillations of the earth of the spheroidal type. The stations are all part of the 
world-wide network of long-period seismographs initially placed in operation during 
the International Geophysical Year and operated by the Lamoiit Geological 
Observatory in cooperation with local institutions. The stations were located at the 
time of the Chilean shock at Agra, India; Hallett, Antarctica; Hong Kong; Lwiro, 
Republic of the Congo; Mt. Tsukuba, Japan; Resolute, Cornwallis Island, North- 
west Territories, Canada; Suva, Fiji; and Uppsala, Sweden. The vertical component 
seismographs all consist of 15-second pendulums and 75-second galvanometers. 

The spectra of these seismograms are presented in this article together with a 
tabulation of the periods of the free oscillations observed at each station. For the 
convenience of the reader all other known data for spheroidal oscillations are also 
tabulated and an average value calculated for each free period. On the basis of these 
data it is possible to compute a phase velocity vs period curve for the fundamental 
Rayleigh mode in the period range of 200 to 3200 seconds. Over much of the period 
range the phase velocities thus obtained have a standard error of the mean of less 
than 0.1 0 . The observed spectral amplitudes are also discussed. 


Discussion of Spectra 

Table 1 lists the sources of the data used to obtain the spectra together with the 
epicentral distance from the stations to the epicenter at 39.5°S, 74.5°W, and also 
the azimuths of the great circle paths between the epicenter and the stations both 
at the epicenter and at the station. The azimuths are given in degrees clockwise 
from 0° at true north. As the origin time of the earthquake was 19:11 :17, 22 May 
1960, table 1 shows that the intervals used for the analysis began between approxi- 
mately four hours to fifteen hours after the earthquake occurred. In most cases the 
galvanometers were driven off scale and were reset the following day when the 
records were (‘hanged. This is the primary reason for the different beginning times 
of the intervals. As much as possible of one complete day’s record was used for 
the analysis. 
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The seismograms were digitized at Lamont using an instrument designed by 
Pomeroy described previously by Alsop, Sutton, and Ewing (1961). All of the 
records were digitized with an interval of approximately six seconds between data 
points. The exact value of the digitization interval was calculated to three decimal 
places by dividing the time length of the record by the number of data points. In 
this way allowance could be made for the stretching or shrinking of the record paper 
during processing. 

Prior to digitizing, the traces were smoothed by eye to remove short period waves, 
i.e., of the order of 10 to 15 seconds. After digitization, in order to reduce the amount 
of data to be analyzed, while avoiding aliasing, the data were smoothed by averaging 
over 19 and 21 points, and then every third point was selected. The data were then 

Table 1 

Summary of Data 


Station 

Latitude 

Longitude 

Epicentral 

Distance 

Azimuth 

from 

Epicenter 

Azimuth 

at 

Station 


Record 


Agra 

27°08'N 

78°01'E 

154. 04° 

110.03° 

234. 64° 

07:31:24, 

23 May - 02:03:38, 

24 May 1960 

Hallett 

72°18’S 

170°18’E 

59. 84° 

198. 66° 

125.93° 

01:42:34, 

23 May - 22:22:45, 

23 May I960 

Hong Kong 22°18'N 

n4°lO'E 

161.32° 

205. 85° 

158.63° 

02:19:00, 

23 May - 01:26:00, 

24 May I960 

Lwiro 

2°15‘S 

28°48'E 

98. 80° 

100. 24° 

229. 63° 

10:40:00, 

23 May - 05:10:24, 

24 May I960 

Mt. Tsukuba 36°1Z*N 

l40°06 f E 

152. 58° 

274. 12° 

107. 42° 

04:44:30, 

23 May - 04:05:3 8, 

24 May I960 

Resolute 

74°4l’N 

94°54*W 

114.70° 

354. 14° 

162.73° 

01:44:41, 

23 May - 01:54:49, 

24 May I960 

Suva 

18°01'S 

178°26’E 

91. 13° 

245. 38° 

132.28° 

23:35:40, 

22 May - 14:33:16, 

23 May I960 

Uppsala 

59°51*N 

I7°38'E 

124.15° 

37. 54° 

249. 12° 

07:46:49, 

23 May - 03:01:40, 

24 May I960 


Fourier analyzed and the periods of the spectral peaks obtained by interpolation. 
This is the same procedure described in slightly more detail by Alsop, Sutton, and 
Ewing (1961). 

The spectra for the eight stations are shown in figures 1 through 8 in increasing 
order of starting time of the seismograms after the earthquake. The periods of the 
observed free vibrations for the fundamental mode are given in the first eight 
columns of table 2. The general appearance of each of the spectra is approximately 
the same. The rapid fall-off of the seismometer sensitivity at long periods limits the 
longest observed free oscillations to those with periods of about 600 to 700 seconds, 
i.e., approximately orders 8 or 9 of the spheroidal mode. The shortest period of 
oscillation observed is approximately 200 seconds. As might be expected, since 
short period oscillations decay more rapidly with time than the long period oscilla- 
tions, the spectra of the seismograms starting many hours after the earthquake 
do not have peaks with as short a period as the spectra of seismograms starting a 
few hours after the earthquake. The general trend of the decay of the short period 
energy may be observed by comparing the spectra in figures 1 through 8 with the 
starting times given in table 1 . The highest spectral peaks occur between 300 and 
400 seconds for the seismograms with an early beginning while for the seismograms 
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Fig. 2. Spectrum of Hallett., Antarctica, seismogram. 
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Fig. 6. Sped nun of Agra, India, seismogram. 
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with a later beginning the highest peaks occur between 400 and 500 seconds or even 
more. The variation of amplitudes between successive peaks will be discussed in 
more detail in another section of this article. 

Spectral peaks corresponding to some of the first higher mode oscillations are 
observed in the spectra for Suva, Mt. Tsukuba, Lwiro, and Uppsala. The oscillation 
iS 8 is observed in all four spectra and there is a suggestion of a peak corresponding 
to this oscillation in the Agra spectrum. (N.B. throughout this article the notation 
of MacDonald and Ness (1961) will be used to identify the various types of free 
vibrations.) The observed first higher mode oscillations are tabulated in table 3. 
The highest order observed is possibly 9. This is in accord with theoretical calcula- 



Fi«. 8. Spectrum of Lwiro, Republic of the Congo, seismogram. 

tions of Alsop (1963) which indicate that the amplitudes of orders 10 through 15 
will be small at the surface of the earth. Theoretical periods for the Gutenberg - 
Bullen B and Gutenberg-Bullen A models from the same article are also included 
in table 3 for comparison. 


Observed Periods 

Table 2 is a composite of all spheroidal free oscillation data published to date for 
the fundamental mode. For orders 10 and below of the fundamental mode there 
can be no confusion in distinguishing spheroidal from torsional oscillations with 
the exception of oS 7 and 0 T 7 so that all of the data available were used for these 
orders. For orders higher than 10 it becomes difficult to positively identify a. spectral 
peak as belonging to a spheroidal or to a torsional oscillation on the basis of period 
alone. Accordingly, for orders higher than 10, only data obtained by vertical seismo- 
graphs or by gravimeters, both of which can only record spheroidal oscillations, 
were used. 

The data selected in accordance with the above criterion, in addition to those 
reported in this article, were taken from the following: Ness, Harrison, and Slichter 
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(1961), designated as Los Angeles in the table; Benioff, Press, and Smith (1961), 
designated as Pasadena; Bogert (1961) designated as Chester; Alsop, Sutton, and 
Ewing (19(31) designated as Palisades although only the data for orders greater 
than seven were obtained on a vertical seismograph at Palisades, the remainder 
being recorded on the first Lamont strain seismograph at Ogdensburg, Xew Jersey; 
Buchheim and Smith (1961) designated as Tiefenort; Bolt (1963) designated as 
Trieste; Connes, Blum, Jobert, and Jobert (1962) designated as Paris; and 
Takeuchi, Saito, Kobayashi, and Xakagawa (1962) designated as Kyoto. All of the 
periods are given in seconds for the purposes of table 2, although originally much of 
the data were given in minutes. 

The average value of these data are also given in table 2. The data contained in 
parentheses were not used in obtaining the average because of the unsatisfactory 
appearance of the peaks in the spectra. The standard error of the mean is also given 


Table 3 

Observed Free Periods of Spheroidal Oscillations, First Higher Mode 



Suva 

Mt. Tsukuba 

Lwiro 

l S 4 



(840.3) 

l s 5 



738.9 

l s 6 



663.8 

1 £ 7 

611.4 


608.8 

1 S 8 

555.7 

553.2 

559.9 

1 S 9 



(514. 6) 


Uppsala 

Average 

Gutenberg 
Bullen B 

Gutenberg 
Bullen A 


(840. 3) 

864 

847 


738.9 

740 

723 


663.8 

666 

649 


‘ 610. 1 

611 

596 

553.2 

555.4 

560 

547 


(514.6) 

511 

501 


for each average value for orders 2 through 37. Above order 37 the data are too 
meager to attempt any statistical analysis. Over much of the range of periods given, 
the standard error of the mean is approximately one part in a thousand. 

The values of period obtained from the analysis of the Resolute seismogram are 
not included in the average. After the averages and standard errors were computed 
it was found that the Resolute values were consistently high. The value of x was 
calculated (Jeffreys, 1948) and was found to be 60.589 for 22 degrees of freedom. 
Tables of x (Eisher and Yates, 1943) indicate that the probability of x being 
greater than 48.268 is less than .001. On this basis the data were rejected as con- 
taining a large systematic error, presently unknown, and the averages and standard 
errors were recalculated. These are the values which appear in table 2. 

Phase Velocities 

Prom the periods of the free spheroidal oscillations given in table 2, the phase 
velocities, 0, of the fundamental Rayleigh mode may be calculated by use of the 
following relation: 


in + \)T 


( 1 ) 
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where a is the radius of the earth, n the order number of the free ose illation, and T 
the period. This is an approximate formula. As discussed in Alsop (1968), for orders 
less than 10 the value of phase velocity obtained is the instantaneous phase velocity 
at an epicentral distance of 90°. 

The phase velocities calculated from the averages of the observed free oscillations 
are listed in table 4 together with the error calculated from the standard error of 


Table 4 

Phase Velocities of Fundamental Rayleigh Mode 



Period 

(sec) 

Phase Velocity 
(km/sec) 


Period 

(sec) 

Phase Velocity 
(km/ sec) 

0 S 2 

3229 

4. 959+ .018 

0 S 23 

315. 5 

5. 399 +. 003 

0 S 3 

2143 

5.337 + . 018 

0 S 24 

306, 1 

5. 338 + . 004 

0^4 

1545 

5.754 + . 011 

0 S 25 

297.6 

5.275 +. 004 

0 S 5 

1189 

6. 121 + .010 

OS 26 

289. 9 

5.211 + . 005 

0 s b 

962 

6. 402 + . 006 

0 S 27 

281.8 

5. 165 + . 006 

0 S 7 

809 

6. 597 + . 008 

0 S 28 

275. 2 

5. 104 + . 006 

0 S 8 

707 

6.661 + . 009 

0 S 29 

268.3 

5.058 + . 006 

0 S 9 

634 

6.646 + . Oil 

0 S 30 

261.9 

5. Oil + .004 

0 S 10 

579 

6. 584 + . 011 

0 S 31 

256. 0 

4. 964 + . 004 

0 S 11 

536 

6.494 + .012 

0 S 32 

250. 1 

4. 925 + . 006 

0 S 12 

502. 3 

6. 375 + . 004 

0 S 33 

245. 1 

4.875 + . 004 

0 S 13 

473. 4 

6.264 + . 004 

0 S 34 

239. 6 

4. 843 + . 004 

0 S 14 

448. 3 

6. 158 + . 004 

0 S 35 

234. 9 

4. 800 + . 006 

0 S 15 

426.5 

6.055 + . 005 

0 S 36 

229. 9 

4. 770 008 

0 S 16 

407. 0 

5. 961 + . 006 

0 S 37 

225. 1 

4. 742 + . 010 

0 S 17 

389. 3 

5.876 + . 005 

0 S 38 

219. 9 

4. 728 

0 S 18 

374. 1 

5.784 +. 005 

0 S 39 

215. 9 

4. 694 

0 S 19 

360. 7 

5.691 + . 005 

0 S 40 

211. 5 

4.673 

0 S 20 

347.4 

5. 621 + .005 

0 S 41 

207. 3 

4.653 

0S21 

335.9 

5. 543 + . 003 

0 S 42 

204. 3 

4.610 

0S22 

325.2 

5.471 + . 005 





the mean of the free period averages. Figure 9 is a plot of the phase velocities as a 
function of period. The phase velocities are obtained of course only at the periods 
of the free oscillations. The curve is achieved by interpolating graphically. At each 
of the points obtained from the free oscillations the error is also indicated. For 
purposes of comparison phase velocity curves obtained from the calculated free 
periods of oscillation for two earth models are also included (Alsop, 1963). The 
models both have velocities according to Gutenberg but one has the densities of 
Bullen’s Model A while the other has the densities of Bullen’s Model B. 

It can be seen from figure 9 that for periods greater than 1 100 seconds the phase 
velocities calculated from the observed free periods of oscillation agree very well 


PHASE VELOCITY (km/sec) 
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with those calculated from the Gutenberg-Bullen B model. For periods less than 
1100 seconds the observed phase velocity diverges from that of the Gutenberg- 
Bullen B model and tends towards the phase velocity curve of the Gutenberg- 



Bullen A model, finally meeting with the curve of the Bullen A model at about 550 
seconds. The observed curve follows the theoretical curve for the Gutenberg-Bullen 
A model, at first slightly below, then slightly above until about 400 seconds. In this 
period range the variations from the theoretical curve are only slight but they are 
consistent and do not show any random character. Below 400 seconds the observed 
phase velocities diverge downwards from the theoretical phase velocities. Once 
again they approach the Gutenberg-Bullen B model but cross over this phase 
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velocity curve and near 230 seconds the observed phase velocities are lower than 
those calculated for either model The observed phase velocities of table 4 in the 
period range of 220 to 340 seconds are in excellent agreement with values obtained 
from surface waves by Brune, Benioff, and Ewing (1901). 

It would appear from the fact that both density models yield phase velocities 
which are higher than the observed phase velocities in the period range of 200-240 
seconds that the fault lies with the velocity model used. This is not too surprising. 
Approximately three-quarters of the path for a free oscillation is oceanic, while 
the Gutenberg velocity model is continental. The “path” of a free oscillation is 
taken to be the path of the two opposing surface waves which form the standing 
wave pattern. 

Using a modification of the program SPHOS 3 described in Alsop (19G3), which 
makes possible the calculation of the free periods of oscillation for earth models 
with a surficial liquid layer, phase velocities have been calculated for the oceanic 
model Case 122 of Sykes, Landisman, and Sat6 (1962). The results which are sum- 
marized in table 5 are disappointing. The phase velocities obtained are even higher 
than those calculated for the two continental models. This model gave excellent 
results for both oceanic Rayleigh and Love waves in the period range of 30 seconds 
to 150 seconds. 

Finally it is of interest to compare the few observed free periods of the first higher 
mode with the theoretical values as shown in table 3. All of the observed free periods 
lie closer to the calculated periods of the Bullen B model than to those of the Bullen 
A model. This is to be expected. Orders 6 through 9 of the fundamental mode do 
have phase velocities more like those of the Bullen A model as the period decreases, 
but the higher mode oscillations have appreciable amplitudes at greater depths 
than those of the fundamental mode. Therefore the higher mode oscillations can 
be expected to favor the same model as do the lower orders of the fundamental mode. 

Spectral Amplitudes 

The vertical amplitudes of the free vibrations of a sphere as a function of position 
on the surface of the sphere vary as P„ m (cos d)e im4> (Morse and Feshbach, 1953), 
where 6 is the colatitude and </> the azimuth angle. For the free vibrations of the 
earth excited by an earthquake, the logical place to choose the pole is at the epi- 
center of the earthquake. For point sources of a simple type — point compressional, 
couple, etc. — it seems reasonable to expect with this choice for the pole that the 
motion for each order n will be characterized by the same degree in or at least by a 
sum of even or of odd m’s of the same symmetry. 

The relative amplitudes of the Legendre polynomials P n were calculated for the 
different epicentral distances given in table 1 using the formula for the asymptotic 
behavior of the Legendre function (Jahnke and Emde, 1945): 

P n(cos 6) = a/ \(l - } } sin 0 — * cot 9 cos 0 | (2) 

y 7T n sin 9 L\ 4 n; 8 n J 

where 0 = (n + £) 6 + tt/ 4. A comparison of the calculated values with the spec- 
tral peaks of the same order showed no similarity. In particular, peaks were found 
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to be missing in the spectra where the Legendre polynomial of that order should 
have appreciable amplitude and were found to be present where the Legendre 
polynomial of that order would have an almost null amplitude. Comparing the 

Table 5 

Free Periods and Phase Velocities 
of Spheroidal Oscillations, Fundamental Mode 
Case 122 



Period (sec) 

Phase Velocity 
(km/ sec) 

00 

£ 

370.9 

5. 832 

0 S 19 

3S6.8 

5.753 

0 S 20 

343.9 

5.678 

0 S 21 

332.1 

5.606 

0 S 22 

321.2 

5. 538 

0 S 23 

311.3 

5. 472 

0 S 24 

302.0 

5.410 

0 S 25 

293.4 

5.350 

0 S 26 

285.3 

5.293 

0 S 27 

277.8 

5.240 

0®28 

270.7 

5.188 

0 S 29 

264.0 

5.139 

0 S 30 

257.7 

5. 093 

0 S 31 

251.7 

5. 049 

0 S 32 

246.0 

5. 007 

0 S 33 

240.5 

4. 967 

0 S 34 

235.4 

4. 929 

0 S 35 

230.4 

4.892 

0 S 36 

225.7 

4. 859 

0S 37 

221.2 

4. 825 

0 S 3S 

216.8 

4.795 

0 S 39 

212.6 

4. 765 

0 

CO 

4* 

O 

208.6 

4.737 

OS4I 

204.8 

4.710 

0 S 4* 

201. 1 

4.684 


missing peaks of the spectra with the predicted zero points of the Legendre poly- 
nomials or the associated Legendre polynomials at the epicentral distance is the 
simplest test which will definitely show if the spectral peaks do not have amplitudes 
similar to polynomials of a given degree m. Therefore the zero point curves for 
Pn and P n 2 of Sato and Usami (1963), who proposed this method, were used to 
determine if the spectra might obey these types of distribution. Again no corre- 
spondence was observed. 
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The magnification of the seismograph as a function of period will cause the ratio 
of the amplitudes to be different from that predicted by considering the relative 
amplitudes of the Legendre polynomials only, but since the magnification curve 
is a smooth curve over these period ranges, it will not introduce peaks and valleys 
into the spectra. Also the time variation of the source must be considered, but for 
periods as long as the free vibrations the spectral characteristics of the source due 
to the time variation should be represented by a smooth curve, most likely 
monotonically decreasing with increasing period. 

In order to explain the difference in phase shift between the vertical and hori- 
zontal motion for different orders of oscillation, Benioff, Press, and Smith (1961) 
postulated a moving source, i.e., faulting with a finite rupture velocity. They found 
that the observed phase shifts were best explained by rupture velocities of 3 to 
4 km/sec and fault lengths of 960 to 1280 km. These authors derived the following 
expression for the amplitude (i n of the vertical motion of the nth order oscillation 
with angular frequency a n due to a moving source: 





PJcos ( 0 


e cos A)] sin 



de 


(3) 


where 0u is the angular fault length, e is the angular variable along the fault, 6 is 
the angular epicentral distance, Co is the rupture velocity, a is the earth’s radius, 
and A is the angle between the great circle through the fault and the great circle 
from the station to the epicenter. 

L T sing expression (3), the relative amplitudes to be expected at the epicentral 
distances given in table 1 were calculated for fault lengths of 960 and 1200 km and 
a rupture velocity of 4 km/sec. Comparison with the spectra showed some similarity 
at the long periods, i.e., greater than 400 seconds, but not at the short periods. 

The discussion up to now has concerned all of the spectra shown in figures 1 
through 8, considered as a group. However, the spectra of Suva, figure 1, is of 
interest by itself. The distance from the epicenter to Suva is slightly more than 90°. 
From (2) it can be seen that near 90° the Legendre polynomials of odd order are 
small for small ft. At the epicentral distance of Suva, they become as large as the 
even order polynomials for n equal to about 37. Formula (2) may be generalized 
for the associated Legendre polynomials (Hobson, 1931) to 


PnX COS 6) 


= (-nr(i 


m(m + 1) 
2 n 


J 2_|Yi + ^_A 

y mr sin 6 L\ 2ft 4ft/ 


•sin (e + t) “ ^ ( m2 “ i) cot 9 cos ( 0 + ?)] 


(4) 


where again © = (ft + %)0 + ir/4. Formula (4) shows that the odd order poly- 
nomials will also be small near 90° for associated Legendre polynomials of even 
degree m, while the even orders will be small for the polynomials of odd degree. 
Figure 1 shows that the odd peaks of the spectra for Suva are indeed smaller than 
the even peaks up to order 21. This tends to indicate that the source has primarily 
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Table 6 

Comparison of Observed Spectral Amplitudes of Suva Spectrum 
with Calculated Amplitudes for Moving Source 



Observed 

Fault Length = 960 km 

Fault Length 

0 S 8 

1.6 

6.4 

8. 8 

0 S 9 

1.2 

1.1 

2.3 

0 S 10 

2.8 

5.9 

7.4 

o s n 

2.2 

1. 3 

2.7 

0 S 12 

2.9 

5. 1 

6.0 

0 S 13 

1.6 

1.5 

3.1 

0 S 14 

4.9 

4.5 

4.9 

0 S 15 

1.2 

1.8 

3.3 

0 S 16 

4.0 

4.0 

4.0 

0 S 17 

1.8 

2.0 

3.5 

0 S 18 

3.6 

3.5 

3.2 

0 S 19 

2.3 

2.1 

3.5 

0 S 20 

3.7 

3.0 

2.7 

0 S 21 

3.8 

2.3 

3.5 

0 S 22 

3.8 

2.6 

2.3 

0 S 23 

1.9 

2.4 

3.4 

0 S 24 

1.9 

2.2 

2. 1 

0 S 25 

2.6 

2.5 

3.2 

0 S 26 

2.7 

1.9 

1.9 

0 S 27 

2. 1 

2.5 

3.0 

0 S 28 

1.8 

1.7 

1.8 

0 S 29 

2.7 

2.5 

2.7 

0 S 30 

1.8 

1.5 

1.7 

0 S 31 

1.8 

2.5 

2.3 

0 S 32 

1.1 

1.4 

1.6 

0 S 33 

1.9 

2.4 

2.0 

0^4 

1. 1 

1.3 

1.4 

0 S 35 

1.2 

2.3 

1.6 

0 S 36 

1.0 

1.2 

1.2 

0 S 37 

0.8 

2.2 

1.2 

0 S 38 

- 

1. 1 

0.9 

0 S 39 

0.9 

2.0 

0.8 


an even symmetry with respect to <t>, the azimuthal angle, since the surface solutions 
are of the form P n m ( cos 6)e ,m *, where m must be the same in both terms because it 
is a separation constant arising in the separation of the scalar wave equation in 
spherical coordinates. The calculations for the moving source show that the effect 
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is to reduce the order at which the amplitudes of the even and odd order poly- 
nomials become equal. For a source length of 960 km and a rupture velocity of 4 
km/sec the even and odd orders should have equal orders at Suva at approximately 
order 23 and for a fault length of 1200 km with the same rupture velocity at ap- 
proximately order 17. Therefore, to this extent, the spectra of Suva, where the even 
and odd order amplitudes become equal at about order 21, could be explained by 
motion along a fault of the dimensions and with a velocity of that suggested by 
Benioff, Press, and Smith. However, further details of the spectrum beyond order 
21 are not in good agreement with the amplitudes predicted for the moving source. 
For purposes of comparison the relative amplitudes predicted for the two different 
fault lengths, for m = 0, are given in table 6 together with the observed amplitudes. 
The calculated amplitudes have been multiplied by a proportionality constant, so 
that the amplitudes for 0 Si 6 are equal. Also the asymptotic value formula for the 
P n 9 s given in (2) has been multiplied by y/n + h to give a normalization in which 
f[P n (x)f dx = 1, where x = cos 8 , (Sommerfeld, 1958). Table 6 shows that for a 
fault length somewhere between 960 and 1200 km the agreement with the observed 
amplitudes is quite good up to order 21, which has a period of 336 seconds, allowing 
for the fact that the sensitivity of the seismograph is low at the long period end of 
the spectrum and consequently the signal-to -noise ratio of the spectrum is small 
for the long periods. There is a trough at n = 23 and 24 which does not occur for the 
calculated amplitudes. Both calculated amplitudes for n = 26 are small; the ob- 
served amplitude is large. Other discrepancies may be noted for higher orders. 

One possible reason for the lack of agreement between the observed and calcu- 
lated amplitudes at the shorter periods for this spectrum and also for the other 
seismograms may be contamination of the spectrum at the shorter periods by some 
of the aftershocks. Another possibility is that while the concept of the moving 
source may be correct, the assumption implicit in (3) that the intensity of the 
radiated energy is uniform all along the fault is open to question. A varying intensity 
of radiation along the fault should have the greatest effect on the short period end 
of the spectrum where the wave lengths are short. 

Conclusions 

Spectral peaks corresponding to the free periods of oscillation of the earth have 
been found in the spectra of eight seismograms obtained from vertical seismographs 
located at various stations around the world. From these and previously reported 
values of the free periods of spheroidal type oscillations, it has been possible to 
obtain very precise values of Rayleigh wave phase velocities in the period range of 
200 to 3200 seconds. While the phase velocity curve agrees with that predicted for 
several earth models over different parts of its period range, it does not agree with 
any existing model over its entire extent. This phase velocity curve should be a 
valuable aid to the problem of modifying existing earth models. Similar data for the 
torsional oscillations would also be valuable and it is hoped that it will be possible 
to obtain these from the horizontal seismograms of these stations. The difficulty 
lies in positively identifying a peak in the spectrum of a horizontal seismogram as 
being definitely a torsional or a spheroidal free period. 

The observed amplitudes of the spectra do not agree well with those predicted 
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for simple point sources. The assumption of a moving source tends to explain the 
observed amplitudes at periods greater than 400 seconds, but not for shorter periods 
than this. It is suggested that variations in the intensity of the radiation emitted 
by the source as a function of time might explain the discrepancies at short periods. 

The peaks of the spectrum of Suva, which are alternately large and small, as is 
to be expected for a station slightly greater than 90° distant from the epicenter, 
indicate that the vertical motion is primarily symmetric with respect to reflections 
through the epicenter. Presumably, the source mechanism must be of this same 
symmetry. 
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